Weakly nonassociative algebras, Riccati and KP hierarchies 
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Abstract 



It has recently been observed that certain nonassociative algebras (called 'weakly nonassociative' , 
WNA) determine, via a universal hierarchy of ordinary differential equations, solutions of the KP 
hierarchy with dependent variable in an associative subalgebra (the middle nucleus). We recall central 
results and consider a class of WNA algebras for which the hierarchy of ODEs reduces to a matrix 
Riccati hierarchy, which can be easily solved. The resulting solutions of a matrix KP hierarchy 
\^ | determine, under a 'rank one condition', solutions of the scalar KP hierarchy. We extend these 

results to the discrete KP hierarchy. Moreover, we build a bridge from the WNA framework to the 
Gelfand-Dickey formulation of the KP hierarchy. 



1 Introduction 



The Kadomtsev-Petviashvili (KP) equation is an extension of the famous Korteweg-deVries (KdV) equa- 
tion to 2+1 dimensions. It first appeared in a stability analysis of KdV solitons [1,2]. In particular, it 
describes nonlinear fluid surface waves in a certain approximation and explains to some extent the forma- 
tion of network patterns formed by line wave segments on a water surface [2]. It is 'integrable' in several 
respects, in particular in the sense of the inverse scattering method. Various remarkable properties have 
been discovered that allow to access (subsets of) its solutions in different ways, see in particular [3-5]. 
Apart from its direct relevance in physics, the KP equation and its hierarchy (see [5, 6], for example) 
is deeply related to the theory of Riemann surfaces (Riemann-Schottky problem, see [7] for a review). 
Some time ago, this stimulated discussions concerning the role of KP in string theory (see [8-11], for 
example). Later the Gelfand-Dickey hierarchies, of which the KdV hierarchy is the simplest and which 
are reductions of the KP hierarchy, made their appearance in matrix models, first in a model of two- 
dimensional quantum gravity (see [12, 13] and references therein). This led to important developments 
in algebraic geometry (see [14], for example). Of course, what we mentioned here by far does not ex- 
^ . haust what is known about KP and there is probably even much more in the world of mathematics and 

physics linked to the KP equation and its descendants that still waits to be uncovered. 

In fact, an apparently completely different appearance of the KP hierarchy has been observed in 
[15]. On a freely generated 'weakly nonassociative' (WNA) algebra (see section [2]) there is a family 
of commuting derivation^] that satisfy identities which are in correspondence with the equations of the 
KP hierarchy (with dependent variable in a noncommutative associative subalgebra). As a consequence, 
there is a hierarchy of ordinary differential equations (ODEs) on this WNA algebra that implies the KP 
hierarchy. More generally, this holds for any WNA algebra. In this way WNA algebras determine classes 
of solutions of the KP hierarchy. 



'Families of commuting derivations on certain algebras also appeared in [16, 17], for example. In fact, the ideas underlying 
the work in [15] grew out of our work in [18] which has some algebraic overlap with [16]. 
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In section [2] we recall central results of [15] and present a new result in proposition [TJ Section [3] 
applies the WNA approach to derive a matrix RiccatH hierarchy, the solutions of which are solutions of 
the corresponding matrix KP hierarchy (which under certain conditions determines solutions of the scalar 
KP hierarchy). In section|4]we extend these results to the discrete KP hierarchy [35-39]. Furthermore, in 
section[5]we show how the Gelfand-Dickey formulation [5] of the KP hierarchy (with dependent variable 
in any associative algebra) emerges in the WNA framework. Section [6] contains some conclusions. 

2 Nonassociativity and KP 

In [15] we called an algebra (A, o) (over a commutative ring) weakly nonassociative (WNA) if 

(a, boc, d) = Va,b,c,d&A, (2.1) 

where (a, b,c) := (aoi)oc-ao(f)oc) is the associator in A. The middle nucleus of A (see e.g. [40]), 

A' := {b G A | (a,b,c) = Va,c G A}, (2.2) 

is an associative subalgebra and a two-sided ideal. We fix / G A, / ^ A', and define a o 1 b := a o b, 

a o n+1 b := a o (/ o n b) - (a o f) o n b , n = 1, 2, . . . . (2.3) 

As a consequence of (12.11 ). these products only depend on the equivalence class [/] of / in A/ A'. The 
subalgebra A(/), generated by / in the WNA algebra A, is called 5-compatible if, for each n G N, 

5 n (f)--=f°nf (2-4) 
extends to a derivation of A(/). In the following we recall some results from [15]. 

Theorem 1 Let A(/) be 5-compatible. The derivations 5 n commute on A(/) and satisfy identities that 
are in correspondence via 5 n i— > dt n ( the partial derivative operator with respect to a variable t n ) with 
the equations of the potential Kadomtsev-Petviashvili (pKP) hierarchy with dependent variable in A'. □ 

This is a central observation in [15] with the following immediate consequence. 

Theorem 2 Let A be any WNA algebra over the ring of complex functions of independent variables 
t\,tz,.... If f G A solves the hierarchy ofODEE 

ft n ■= d tn (f) = f °n / , n = l,2,..., (2.5) 

then —fa lies in A' and solves the KP hierarchy with dependent variable in A'. □ 

Corollary 1 If there is a constant v G A, v A', with [v] = [f] G A/A', then, under the assumptions 
of theorem\2\ 

4> := v - f E A' (2.6) 

2 Besides their appearance in control and systems theory, matrix Riccati equations (see [19-21], for example) frequently 
showed up in the context of integrable systems, see in particular [22-34]. 

3 / has to be differentiable, of course, which requires a corresponding (e.g. Banach space) structure on A. The flows given 
by J2.5b indeed commute [15]. Furthermore, i2.5l implies 5-compatibility of the algebra A(/) generated by / in A over C [15]. 
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solves the potential KP (pKP) hierarch^ 

3 

i j',fc=l 

where is totally antisymmetric with £123 = 1, Aj, ? = 1,2,3, are indeterminates, and <f>±\\](t) := 
<j>{t ± [A]), where t = (ti,t 2 , . . .) and [A] := (A, A 2 /2, A 3 /3, . . .)• □ 

Remark 1. If C € A' is constant, then / = v' — (<ft + C) with constant v' := v + C satisfying 
[z/] = [u] = [/]. Hence, with also + C is a solution of the pKP hierarchy. This can also be checked 
directly using (12.71 ). of course. □ 

The next result will be used in section [4] 
Proposition 1 Suppose f and f solve ( 12.51 ) and [f] = [f] in a WNA algebra A. The equation 

f'o f = a (f-f) (2.8) 

is then preserved for all a S C. 

Pm0f: (/'°/k = fl n of + fof tn = (fo n f')of + f'o(fo n f) 
= (/' °n /') o / - /' o n (/' o /) + a f o n (/' - /) 
+/'<>(/ °n /)-(/' of) o n f + a (f-f)o n f 
= ~f °n+i / + /' o n+ i / + a {f' o n f - f o n f) = a (/' - f) tn ■ 

In the third step we have added terms that vanish as a consequence of (12- 8b - Then we used (3. 1 1) in [15] 
(together with the fact that the products o n only depend on the equivalence class [/] = [/'] S A/A'), and 



also (12- 3b . to combine pairs of terms into products of one degree higher. □ 
Remark 2. In functional form, (12.51) can be expressed (e.g. with the help of results in [15]) as 

A" 1 (/-/- [ A])-/-[A]°/ = 0. (2.9) 
Setting /' = /_[a] (which also solves (12.91 ) if / solves it), this takes the form (12.81) with a = —A" 1 . □ 



In order to apply the above results, we need examples of WNA algebras. For our puiposes, it is 
sufficient to recall from [15] that any WNA algebra with dim (A/ A') = 1 is isomorphic to one determined 
by the following data: 

(1) an associative algebra A (e.g. any matrix algebra) 

(2) a fixed element g € A 

(3) linear maps £,1Z : A — > A such that 

[C,K]=0, C(a o b) = C{a) o b , K(a o b) = a o K(b) . (2.10) 

Augmenting A with an element / such that 

f°f'=9, foa:=C(a), aof:=n(a), (2.11) 

leads to a WNA algebra A with A' = A, provided that the following condition holds, 

3a,b £ A : lZ{a) ot/ao C(b) . (2.12) 

This guarantees that the augmented algebra is not associative. Particular examples of C and 1Z are given 
by multiplication from left, respectively right, by fixed elements of A (see also the next section). 

4 This functional representation of the potential KP hierarchy appeared in [41,42]. See also [15,26] for equivalent formulae. 
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3 A class of WNA algebras and a matrix Riccati hierarchy 

Let A4(M, N) be the vector space of complex M x N matrices, depending smoothly on independent 
real variables t\, t%, . . ., and let S, L, R, Q be constant matrices of dimensions M x N, M x M, N x N 
and N x M, respectively. Augmenting with a constant element v and setting^] 

vov = -S, v o A = LA , Aou = -AR, Ao B = AQB , (3.1) 

for all A, B e M(M, N), we obtain a WNA algebra (A, o). The condition (12.121 ) requires 

RQ + QL. (3.2) 

For the products o n , n > 1, we have the following result. 

Proposition 2 

z, o n z/ = -S n , vo n A = L n A, Ao n v = -AR n , Ao n B = AQ n B, (3.3) 



where 



f n Q L n ) =Hn Witk H ' = ( S L 1 (3 - 4) 



Proof: Using the definition (12.31 ). one proves by induction that 

QL n + RQ 

for n = 1, 2, . . ., where Si = S, L\ = L, Ri = R, Q\ = Q. This can be written as 

Rn+l Qn+1 \ I Rn Qn 



which implies (13.4I ). □ 
Using (12.61 ) and (13.31 ) in (12.5b . leads to the matrix Riccati equation^] 

<t>t n = S n + L n (j> - <t>Rn - <t>Q n <f> , 71=1,2,.... (3.5) 

Solutions of (13.51 ) are obtained in a well-known way (see [21,32], for example) via 

(/) = YX- 1 (3.6) 

from the linear system 

Z tn = H n Z, Z =(y) (3 " 7) 
with an N x N matrix X and an M x N matrix Y, provided X is invertible. This system is solved by 
Z(t) = e^ H) Z where £(H) :=^t n H n . (3.8) 

n>l 



5 Using J2.6b . in terms of / this yields relations of the form J2.1U . 

6 The corresponding functional form is A _1 (</!>[ A ] —<£>) + 4>Q4>[\] = S + Lcf>\x\ — 4>R, which is easily seen to imply i2.7l , 
see also [43]. The appendix provides a FORM program [44,45] which independently verifies that any solution of J3.5b . reduced 
to n — 1, 2, 3, indeed solves the matrix pKP equation in (Ai(M, N), o). 
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If Q has rank 1 , then 



<p := HQcj,) (3.9) 

defines a homomorphism from (M.(M, N), o) into the scalars (with the ordinary product of functions). 
Hence, if 4> solves the pKP hierarchy in (Ai(M, N), o), then tp solves the scalar pKP hierarchy^] More 
generally, if Q = VU T with V, U of dimensions N x r, respectively M x r, then U T (j)V solves the 
r x r-matrix KP hierarchy. 

GL(N+M, C) acts on the space of all (N+M) x (N+M) matrices H by similarity transformations. 
In a given orbit this allows to choose for H some 'normal form', for which we can evaluate (13.81 ) and then 
elaborate the effect of GL(N + M, C) transformations (see also remark 3 below) on the corresponding 
solution of the pKP hierarchy, with the respective Q given by the normal form of H. By a similarity 
transformation we can always achieve that Q = and the problem of solving the pKP hierarchy (with 
some non-zero Q) can thus in principle be reduced to solving its linear part. Alternatively, we can always 
achieve that 5 = and the next two examples take this route. 

Example 1. If 5 = 0, we can in general not achieve that also Q = 0. In fact, the matrices 

l) and H :=(« a ») ,3,0, 

are similar (i.e. related by a similarity transformation) if and only if the matrix equation Q = RK — KL 
has an N x M matrix solution K [56-60], and then 

« = ™„t-., r=('--M. <3,i) 



It follows that 



and thus 



fl» = THST-i = ( ^ RnK ~ KLn ) (3.12) 



j {B) = i " n " '"mr ) . (3.i3) 





If (13.21 ) holds, we obtain the following solution of the matrix pKP hierarchy in (Ai(M, N), o), 

$ = e ^ L Uo (In + K<f> - Ke^Uo^e'^ , (3.14) 

where = YqXq . This in turn leads to 

(p = tr 

= tr (log(Ijv + Kcf) - e- m Ke^ L) (j) ) 
= (logr) tl , T:=det(I N + K<p -e-^Ke^cPo) ■ (3.15) 



If rank(Q) = 1, then tp solves the scalar pKP hierarchy. Besides (13.21) and this rank condition, further 
conditions will have to be imposed on the (otherwise arbitrary) matrices R, K, L and <pQ to achieve that ip 



7 For related results and other perspectives on the rank one condition, see [46] and the references cited there. The idea to 
look for (simple) solutions of matrix and more generally operator versions of an 'integrable' equation, and to generate from it 
(complicated) solutions of the scalar equation by use of a suitable map, already appeared in [47] (see also [48-55]). 
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is a real and regular solution. See [61], and references cited there, for classes of solutions obtained from 
an equivalent formula or restrictions of it. This includes multi-solitons and soliton resonances (KP-II), 
and lump solutions (passing to KP-I via t2 n l— * i tin an d performing suitable limits of parameters). 
Example 2. If M = N and S = 0, let us consider 

H = TH T~ l , H =(q T=f I Q ~fY (3.16) 

with I = Jjv and a constant N x N matrix K. As a consequence, 

Q = I+[L,if]. (3.17) 
We note that Hq is not similar to diag(L, L) [60]. Now we obtain 



IT -7)11,1 ' - j L Q " + ) 13.181 



and furthermore 

_ ( e 

e^) 

If [L, [L, if]] 7^ (which is condition (13.21 )). we obtain the solution 

= e ?(L Vo (/ + ^^o + F)- l e~ i{L) (3.20) 
of the matrix pKP hierarchy in (M(N, N), o), where 

F := (^n^L™- 1 - e^ife^)^ ■ (3.21) 

n>l 

Furthermore, using F tl = e~^ L \l + [L, K])e^ i Vo> we find 

^ = tr((/+[L,F])</») =tr(F tl (/ + ^ + F)- 1 ) = (tr log(7 + Kfo + F)) tl (3.22) 

and thus 

99 = (logr) tl , T-detfj + Kfo + i^ntnL^-e-^Ke^Wo) . (3.23) 

n>l 

If rank(i + [L, i£]) = 1 (see also [46,62,63] for appearances of this condition), then <p solves the scalar 
pKP hierarchy. Assuming that <fto is invertible, we can rewrite r as follows, 

r = det (e^ + K)e~^ +^nt n L n_1 - if) (3.24) 

n>l 

(dropping a factor det(<fo))- This simplifies considerably if we set <fi = — Choosing moreover 

Lij = ~(qi - qj)~ l i^j, La = -pi , K = diag(<?i, ...,q N ), (3.25) 

(13.241 ) reproduces a polynomial (in any finite number of the t n ) tau function associated with Calogero- 
Moser systems [46,62,63]. Alternatively, we may choose 

L = diag(<?i, ...,q N ), Kij = (g* - g,) -1 i± j , K u =pi . (3.26) 



8 Note that in this case 4> = (X] n >i n L n 1 — K) 1 , which is rational in any finite number of the variables t n . 
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The corresponding solutions of the KP-I hierarchy fan > itm) include the rational soliton solutions 
('lumps') originally obtained in [64]. In particular, N = 2 and q2 = —qi, P2 = Pi (where the bar means 
complex conjugation), yields the single lump solution given by 



T = l^+^l)| 2 + 4^)2 where £(<!)■= Y,™t n q n - 1 



(3-27) 

{t 2k ^it 2 k ,k=l,2,...\ 



Example 3. Let M = N and L = Stt^, R = tt + S, Q = 7r + S'7r_, with constant N x N matrices 
7r + , 7r_ subject to 7T + + 7T_ = I. The matrix H can then be written as 



which lets us easily calculate 



(3.28) 



H n = ( ^tf " ^tf ) . (3.29) 



As a consequence, we obtain 

= (-C+ + e €(5) C_)(vr_C + + TT+e^C-) -1 , (3.30) 

where C± := I T 7 ^±4'o- This solves the matrix pKP hierarchy in Ai(M, N) with the product Ao B = 
Att+Stt-B if (13.21 ) holds, which is tt+S(tt + — 7r_)S"7r_ ^ 0. If furthermore rank(7r + 5-7r_) = 1, then 

^ = tr(Q0) = -tr(7r + 5) + (logT) tl , r = det(ir_C + +ir + ^C_) (3.31) 

solves the scalar pKP hierarchy (see also [43]). We will meet the basic structure underlying this example 
again in section [5] 

Remark 3. A GL(N + M, C) matrix 

T=(il\ (332) 



C D 



can be decomposed as follows, 



H';r)( s ;:)(A:^ 



if D and its Schur complement Sp = A — BD~ X C are both invertible. Let us see what effect the three 
parts of T induce on when acting on Z. 

(1) Writing P = D~ l C, the first transformation leads to i— > + P, a shift by the constant matrix P. 

(2) The second transformation amounts to ^ ^ D^S^ 1 (where (j> is now the result of the previous 
transformation). 

(3) Writing K = —BD^ 1 , the last transformation is (f> h-> (f> ( I N - Kcj))^ 1 . □ 



4 WNA algebras and solutions of the discrete KP hierarchy 

The potential discrete KP (pDKP) hierarchy in an associative algebra (A, o) can be expressed in func- 
tional form as follows£j| 

n(A) + - n(A)_ M = ^Gu) + - n(M)_[A] , (4-D 



This functional representation of the pDKP hierarchy is equivalent to (3.32) in [39]. 
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where A, \i are indeterminates, 

:= \-\</> - <j>_ [x] ) - (4> + - 0_ [A] ) o 4> , (4.2) 

and = (4>k)kez, </>t '■= 4>k+l- The pDKP hierarchy implies that each component fa, k G Z, satisfies 
the pKP hierarchy and its remaining content is a special pKP Backhand transformation (BT) acting be- 
tween neighbouring sites on the linear lattice labeled by k [35, 39]. This suggests a way to extend the 
method of section [3] to construct exact solutions of the pDKP hierarchy. What is needed is a suitable 
extension of (12.51 ) that accounts for the BT and this is offered by proposition Q] 

Theorem 3 Let Abe a WNA algebra with a constant element v € A, v A'. Any solution 

f = (y - fa) k& , (4-3) 

of the hierarchy ( 12.51 ) together with the compatible constraint 

/+ o / = (4.4) 

yields a solution (ft = (fa)kez of the pDKP hierarchy in A'. 

Proof: Since [/ + ] = [/], the compatibility follows by setting /' = / + and a = in proposition Q] Using 
ft 1 =f° f, we rewrite (12.91 ) as 

A _1 (/ ~ /-[A]) + (/ " /— [A] ) O / - f tl = . 

Inserting / = v — <p, this takes the form 

a _1 o - cf)_ [x] ) - 4> tl -(</>- o4> = e- 6L [A ] 

with 9 := —(f) o v. Next we use (14.41 ) and f tl = f o f to obtain (/+ — /) o / + f tl = 0, which is 

fa, - (0+ - 4>) o $ = 9+ - 9 . 
Together with the previous equation, this leads to 

A- 1 ^ - fa [x] ) ~ - <t>-[\\) o(t> = 9 + - 9_ [x] 
(which is actually equivalent to the last two equations), so that 

n(X) = 9 + -9_ [x] . 



This is easily seen to solve (14. II ). □ 

Let us choose the WNA algebra of section l3f^l Evaluation of (12.51 ) leads to the matrix Riccati 
hierarchy (I3.5I ). and (14.41 ) with f + = v + C — <j) + becomes 

S + CR+(L + CQ)(/)-<p + R-<p + Q(l) = 0, (4.5) 

which can be rewritten as 

<P + = (S + L<t>){R + Qcpy 1 + C = Y + (X+)- 1 (4.6) 



10 Note that @3 implies f n+ o„ / = 0, where /"+ := f k+n . This follows by induction from f<. n +V+ o n+1 f = 

j(n + l)+ Q Qn _ (J(n + 1)+ yn+) ^ j = y(n + l) + Q ^ ^ _ ^+ Q ^ ^ where wg use( j jjjj 

[/ n+ ] = [/] in the nrst step- 

"Since there is only a single element v, the matrices L, R, S do not depend on the discrete variable k. 
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(assuming that the inverse matrices exist), where X + ,Y + are the components of 

Z + = TH Z = TH e^ H) Z^ , (4.7) 

with Z, H, T taken from section [3] Deviating from the notation of section [3j we write Z^ for the 
constant vector, since Zq should now denote the component of Z at the lattice site 0. In order that (14.71 ) 
defines a pDKP solution on the whole lattice, we need H invertible. Since the matrix C, and thus also 
T, may depend on the lattice site k, solutions of (14.11 ) are determined by 

Z k = T k HT k ^H ■ ■ ■ T X HZ Q , Z_ k = (T^H)' 1 (T^H)- 1 ■ ■ ■ {T^H)- l Z Q , k G N . (4.8) 

This corresponds to a sequence of transformations applied to the matrix pKP solution <po determined by 
Zq, which generate new pKP solutions (cf. [35]). (pi is then given by (14.61 ) in terms of <pQ, and 

<fa = [LS + SR + LdR+il? + SQ + LC 1 Q)(j> ] 

x [R 2 + QS + QdR + (QL + RQ + QdQ)^]' 1 + c i ( 4 - 9 ) 

shows that the action of the T k becomes considerably more involved for A; > 1. In the special case 
T k = In+m (so that C k = 0), we have 

Z k = e^ H) {H k zf ] ) k£Z. (4.10) 

If x[ , are the components of the vector H k Z^°\ the lattice component <\> k of the pDKP solution 
determined in this way is therefore just given by the pKP solution of section [3] with initial data (at t = 0) 

0f = Y^\ X f r 1 = L k 4>f [R k + (R k K - KL k )4>Q > Y~ 1 . (4.11) 

With the restrictions of example 1 in section |3j assuming that L and R are invertible (so that H is 
invertible), the corresponding solution of the matrix pDKP hierarchy (in the matrix algebra with product 

A o B = A(RK - KL)B) is 

4> k = e ^L k ^[R k (I N + Kcpf) - e-tWK<*ML k <j>f]- x e-W , k G Z , (4.12) 

which leads to 

<p k = (logTfc)^ with r k = det (R k {I N + K(P [ q } ) - e-^ K) Ke^ L) L k (^ ) \ keZ. (4.13) 



If Q = RK — KL has rank 1 , this is a solution of the scalar pDKP hierarchyf^l As a special case, let us 
choose M = N, L = diag(pi, . . . ,Pn), R = dia,g((/i, . . . , gjv), an d K with entries Kij = p j ) 1 IlH 
Then Q has rank 1 and we obtain iV-soliton tau functions of the scalar discrete KP hierarchy. These can 
also be obtained via the Birkhoff decomposition method using appropriate initial data as in [65, 66]. 

With the assumptions made in example 2 of section [3j setting 4>^ = —K' 1 , assuming that K and 
L are invertible, and choosing for T k the identity, we find the matrix pDKP solution 

<f) k = (S^ntnL 71 - 1 + kL~ l - K) , k£Z. (4.14) 

n>l 

If rank(ijv + [L, K]) = 1, this leads to the following solution of the scalar pDKP hierarchy, 

<Pk = 0ogTk)t 1 with r fc = det(^nt n L n - 1 + /cL- 1 - k} . (4.15) 

n>l 

In example 3 of section |3l H is not invertible, so that (14.71 ) does not determine a pDKP solution. 



12 Recall that tp = tr(Q</>) (cf. 13.9b determines a homomorphism if Q has rank 1. As a consequence, if 4> solves the matrix 
pDKP hierarchy | |4.U . then ip solves the scalar pDKP hierarchy. 



13 The condition d3.2fc requires qt ^ pj for all i,j = l,...,N. 
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5 From WNA to Gelfand-Dickey 

Let 1Z be the complex algebra of pseudo-differential operators [5] 

V=Y i v i &, (5.1) 

i<t.oo 

with coefficients Vi G 21, where 21 is the complex differential algebra of polynomials in (in general 
noncommuting) symbols ul™', m = 0, 1, 2, . . ., n = 2, 3, . . ., where d(u^) = ui m+1 ' and d(vw) = 
d[v) w + v d(w) for v, w G 21. We demand that ti« m ' ) , n = 2, 3, . . ., m = 0, 1, 2, . . ., are algebraically 
independent in 21, and we introduce the following linear operators on 1Z, 

S(V) := £V , 7T+(V) := V> , vr_(V) := V< := V - V> , (5.2) 

where V>o is the projection of a pseudo-differential operator V to its differential operator part, and 

£ = d + u 2 d- 1 + u 3 d- 2 + ■ ■ ■ . (5.3) 

Let / denote the identity of 1Z (which we identify with the identity in 21), and let O be the subspace of 
linear operators on 1Z spanned by S and elements of the form S7r±S7r± ■ ■ ■ ir±S (with any combination 
of signs). O becomes an algebra with the product given by 

AoB:= Att+Stt-B . (5.4) 

(O, o) is then generated by the elements (Sir-) m S (n + S) n , m, n = 0, 1, Let us furthermore intro- 
duce A := {v G 21 : v = res(^4(/)), A G O}, where res takes the residue (the coefficient of d^ 1 ) of a 
pseudo-differential operator. This is a subalgebra of 21, since for A,5eOwe have 

res(A(/)) res(B(J)) = res(^7r +( S^_B(/)) , (5.5) 

so that the product of elements of A is again in A. As a consequence of this relation (read from right to 
left), A is generated by the elements res((S , 7r_) m S'(7r + S') ra (/)), m, n = 0, 1, — Based on the following 
preparations, we will argue that A and (O, o) are actually isomorphic algebras. 

Lemma 1 For all V G 1Z, 

res((,S7r_) m V) = res(P m V) , m = 0, 1, . . . , (5.6) 
where T>q = I and {V m }^ =1 are the differential operators recursively determined by V m = (Z? m _i£)>o. 
Proof: We do the calculation for m = 2. This is easily generalized to arbitrary m G N. 

res((^_) 2 V) = res(£(£V< )<o) = res(£> £V<o) = res((£> £)> V) = res(V 2 V) . n 

Proposition 3 

m / \ 

res((S , 7r_) m S'(7r + S') n (/)) = ( j u^ +n+2 _ k + terms nonlinear in , m, n = 0, 1, . . . (5.7) 

fc=o ^ ^ 
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Proof: According to the preceding lemma, we have 

res((Sir_) m S(ir + ST(I)) = ies(V m S(ir + S) n (I)) . 

Next we note that V m = d m + D m , (ir + S) n (I) = d n + D' n with differential operators D m ,D' n (of 
degree smaller than m, respectively n) such that each of its summands contains factors from {u^ } (so 
their coefficients are non-constant polynomials in the u k , ). It follows that 

res((57r_) m 5(7r + 5)"(/)) = res((5 m + D m )£ <0 (d n + D' n )) 

= res(<9 m £<o<9 n ) + terms nonlinear in u% . 

It remains to evaluate 

oo oo m / \ 

res(5 m £ <0 ^) = £ re S (d m u 1+3 d n ^) = ]T res( ^ I ™ ) ug^ m+n 
j=i i=i fc=o ^ ' 

771 

= £ 



fc=0 



(fc) 

Tj, )' U m+n+2-fc • □ 



According to the last proposition, the linear term with the highest derivativd_j in the residue of 
(Sir-) m S(Tr + S) n (I) is given by uQ 2 . We conclude that the monomials (STT^) m S(-K + S) n , m,n = 
0, 1, . . ., are algebraically independent in (O, o), since any algebraic relation among them would in- 
duce a corresponding algebraic relation in the set of u n , but we assumed the u n to be algebraically 
independent. Together with (15.5I ). this implies that A and (0, o) are isomorphic algebras. 

The last result allows us to introduce a WNA structure directly on A as followslll Augmenting A 
with / such that, for V, W £ 0{I), 

fof:= -res(fi) , / o res(V) := res(£V <0 ) , 

res(V) o / := -res(V< £) , res(V) o res(W) := res(V) res(W) , (5.8) 

indeed defines a WNA algebra A = A(/). The relations (15.81 ) are well-defined since res(vl(/)) uniquely 
determines A € O. By induction we obtain 

fo n f = -res(£ n ) , / o n res(V) = res(£ n V< ) , 

res(V) o n f = -res(V< £ n ) , res(V) o n res(>V) = res(V< £ n >V<o) • (5.9) 

Let the u n now depend on variables i 1; i 2 , . . ., and set <9 = <9 4l . The hierarchy (12.51 ) of ODEs, 

ft n = f o n f = -res(£ n ) , n = l,2,..., (5.10) 

by use of the WNA structure implies 

a tn (res(£ m )) = -d tn (f°mf) = -ft n °mf-fo m f tn 

-(/ on /) ° m f ~ f ° m (/ °n /) = re S (z m (2 n )<o - (£ n )>o£ m ) 

= res([(n> ,£ m ]) • 



14 If m = 0, the linear term is simply u n +2 and thus again 'the linear term with the highest derivative' 
15 Note that the corresponding WNA structure for (0, o) resembles that of example 3 in section[3] 
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Since also d tn {res(£ m )) = res([(£ m )> , £"]) = <9 tm (res(£ n )), we conclude that if we extend A to A 
by adjoining an element <p = d^ 1 {u2), then 

<j) tn = res(£ n ) , n = 1,2,... . (5.12) 

It follows that i> := f + 4> satisfies dt n (y) = 0, n = 1, 2, . . ., and is therefore constant. (15.121 ) determines 
all the Ufc in terms of the derivatives of (j) (see [67], for example). From (15.121 ) with n = 2, 3, and (15.1 II ) 
with m = n = 2, we recover the pKP equation 

(40t 3 - &ititi - S^t! 2 )^ -3<t>t 2 t 2 + 6 [^ti)^] = 0, (5.13) 

in accordance with the general theory. More generally, the equations (15.111 ) determine the whole pKP 
hierarchy. They are the residues of 

d tn (£ m ) = [(£")> ,£ m ], m,n=l,2,.... (5.14) 

This is equivalent to the Gelfand-Dickey (GD) system dt n (£) = [(£ n )>o, £], n = 1, 2, . . ., which is a 
well-known formulation of the KP hierarchy (see [5], for example). 

We have thus shown how the Gelfand-Dickey formulation of the KP hierarchy can be recovered in 
the WNA framework. In fact, for the particular WNA algebra chosen above, the hierarchy (12.51 ) of ODEs 
is equivalent to the Gelfand-Dickey formulation of the KP hierarchy. 

6 Conclusions 

In this work we extended our previous results [15, 61] on the relation between weakly nonassociative 
(WNA) algebras and solutions of KP hierarchies to discrete KP hierarchies. We also provided further 
examples of solutions of matrix KP hierarchies and corresponding solutions of the scalar KP hierarchy. 
In particular we recovered a well-known tau function related to Calogero-Moser systems in this way 
(example 2 in section©. Furthermore, we established a connection with the Gelfand-Dickey formulation 
of the KP hierarchy. As a byproduct, in section|5]we obtained a new realization of the. free WNA algebra 
generated by a single element, which also has a realization in terms of quasi-symmetric functions [15]. 
There is more, however, we have to understand in the WNA framework. In particular this concerns the 
multi-component KP hierarchy (see [68] and references therein) and its reductions, which include the 
Davey-Stewartson, two-dimensional Toda lattice and iV-wave hierarchies. Our hope is that also in these 
cases the WNA approach leads in a quick way to relevant classes of exact solutions. 

Appendix: From Riccati to KP with FORM 

The following FORM program [44,45] verifies that any solution of the first three equations of the Riccati 
hierarchy (13.51 ) solves the pKP equation in an algebra with product A o B = AQB. 

Functions phi, phix, phiy, phit, L, Q, R, S, dx, dy, dt; Symbol n; 

Local pKP = dx*(4*phit - 6*phix*Q*phix - dx~2*phix) - 3*dy*phiy 

+ 6* ( phix*Q*phiy - phiy*Q*phix ) ; * pKP equation 

repeat; 

id phix = S + L*phi - phi*R - phi*Q*phi; * Riccati system 

id phiy = S(2) + L(2)*phi - phi*R(2) - phi*Q ( 2 ) *phi ; 

id phit = S(3) + L(3)*phi - phi*R(3) - phi*Q ( 3 ) *phi ; 

id dx*phi = phix + phi*dx; * product rule of differentiation 
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id dy*phi = phiy + phi*dy; id dt*phi 
id dx? { dx, dy, dt } *L? {L, Q, R, S } = L*dx; 



phit + phi*dt; 

* L,Q,R, S are constant 



* recursion relations for matrices (see proof of proposition 2) : 

id L(n?{2,3}) = L*L(n-l) + S*Q(n-l); 

id R(n?{2,3}) = Q*S(n-l) + R*R(n-l); 

id S(n?{2,3}) = L*S(n-l) + S*R(n-l); 

id Q(n?{2,3}) = Q*L(n-l) + R*Q(n-l); 

id L?{L,Q,R,S} (1) = L; 

endrepeat; 

id dx? { dx, dy, dt } = 0; 

print pKP; * should return zero 



This program provides an elementary and quick way toward the classes of exact solutions of the KP 
equation given in the examples in section [3] 
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